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SUMMARY 


OBJECTIVE 


Determine the high-frequency acoustic scattering characteristics of a random set of ice 
keels. 


APPROACH 


Parabolic sections were used to represent ‘ce keels that are random in size and spacing. 
A large number of rays were compuicd tc model the interaction of a plane wave with the 
random set. 


RESULTS 


The principal results of this work are diagrams of scattering strength versus angle as 
functior.s of incident angle. Reflection lo:s data in the specular direction are also given. 


CONCLUSIONS 


|. Plane wave reflection coefficients strongly control the scattering action of the set 
of keels. 


At grazing angles below 20 degrees, transmission through the keels as shear waves 
is the dominant energy path. 


tw 


3. Specular reflection is determined predominantly by the extent of unshadowed 
surface. 


4. The parabolic-shaped keels underestimate backscatter from the keels. 


RECOMMENDATIONS 


1. Determine and paramaterize the shapes and internal structure of real keels, and if 
substantial variation is observed, determine probability distributions for these 
parameters. 


2. Use compound figures to mode! realistic keel shapes and their acoustic properties. 
3. Include internal structure and attenuation in models. 
4. Seek experimental verification for acoustic results reported here. 


. Aceessien Por 
WTIS GRaal rd 
3 QO 


Distribution/ 
Availability Codes 
lavatl and/or 


Ae 


CONTENTS 


SU NAIR AN RGVar oh ses Noes a dl ot end meet clio a eile earn S, itt 
INTRODUGHIONED te a eee Ck SSRN | Mien Be Fore 
COMPUTER ET OD ee ete ane dt 1 
SOMTIER ING RESUS ee a eee eu 5 
GEOMETRIC IMPLICATIONS OF SCATTERING ANGLES ........-20-c00e000c00e 9 

BACKSCATHERW en cmt nol Liiw oietT_supanitom die nate, ON | 9 

SURFAGHSHADOWING 8 oie eee ee 10 

SCATTERING EIINGTLONS sek! euler, 21taek uabwe les bemmtnde we So 12 
MOD ELGSELORAKG O MTENIGS i whe os esha ays tones oe de cme NG 15 
GONCLUSIONS. eee a. tpi OA IT MR He ROAR ICS NG ote 16 
RECOM MENDATIONS Mtns vii noun ae one nN ay ee eomR Ee 16 


ha oe 


‘ 


LP ae a a) 


~ 


eh, Sa 


at 
a 
# 
< 


alias 98 a 
PP Pe 
rat 


bg ee 


INTRODUCTION 


The scattering of rays by a random collection of ice keels has been investigated by the 
Naval Ocean Systems Center (NOSC). This repert presents the results of this investigation. The 
computer ray program used in this project was reported at the Austin meeting of the Acoustic 
Society of Amenca in Apnil 1985 and in Ref. 1. 


In a 1976 article in the Jcurnal of the Acoustic Society of America (Ref. 2), 
O.I. Diachok showed that Twersky’s boss-scattering theory was an effective approach to 
undenice scattering and propagation. A first-term approximation from that theory was shown 
to give effective rcsults for high frequencies, as did another at low frequencies. The current task 
derives high-frequency numerical values similar to the former, the main difference being that 
the individual bosses are of random size. 


This report wiil first discuss the random ice keel model used, and then the ray 
computation technique. There will then be a discussion of two functions that were derived or 
estimated, the anguiar distribution of rays scattered from the set of keels and the specular 
scattering coefficient The final section wil! discuss some of the geomeiry and mechanics of the 
computations to explain vanous observed features of the computed functions. These discus- 
sions will aid in determining which results arise from given assumptions of the modei and 
provide a basis for estimating the acoustic oroperties of real keels. Such estimations, however, 
are not attempted here. An alternative view of ice keels, in which their rough or ble ck 
character is emphasized, has been provided by S. Chin-Bing (Ref. 3). 


COMPUTER METHOD 


Figure | indicates the basic computationai strategy. A large number of equally spaced 
rays are launched towards a set of keels represented as parabolas. As the rays emerge from the 
set of keels, they are accumulated in angular bins. The contents of these bins then represent the 
scattering cnaractenistics of the keel set. The distance between keels is random, drawn from a 
rectangular distribution to give 9.5 keels; km. The depths of the keels are drawn from a 
Rayleigh distnbution with a mean depth of 4 m. A further parameter, the ratio of depth to 
width is | to 3.2. Therefore, the kee!s are all similar geometric figures. Most of these para- 
meters are taken from Diachok’s 1976 article. 


Figure 2 shows details of an individual keel on an equal length-depth scale. This ts the 
cross section of an elongated object. For simplicity the keel is assumed to be infinite in length. 
The maximum slope of the keel edge is at its upper edge and is 51.3 degrees. The keel model as 
shown. was next modified to simulate an orientation of the keels that is random when viewed 
from above. This was done by increasing the width of the keels by a function of a random 
angie between 0 and 90 degrees — approximately the secant of the angle. However, this 
function is modified by the fact that, in a given range interval, one is more likely to cross a keel 
perpendicular to a path than one nearly parallel. 
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Figure 1. Model of a random sei of ice keels, with incoming and 
scattered rays suggested by the arrows. 


For random onentation 


Mu.tipty wedth by (1-2) 
where x is random. 0 <a <! 


Figure 2 Cross section of a keel model The equation gives the 
widening factor caused by oblique incidence on 3a randomly 
oriemed keel 


If our line of sight is crossed by a keel of width w at an angle of incidence of 6, then the 
apparent width of the keel is sec 6. The probability of crossing a randomly oriented keel at 
angle 6 must he derived. 


This probability is proportional to cos 6. By symmetry, @ can be confined to the 
interval from 0 to 90 degrees. The probability distribution function is therefore cos 0, 09<6<90 
degrees and zero elsewhere. The cumulative distribution function of 6 is therefore sin 6 over the 
same interval. A random value of 6 from this distribution can now be picked by taking arcsin 
(x) where x is uniformly random from 0 to |. The desired width of this randomly selected keel, 
W, is then w sec [arcsin(x)] or 


W=w(l - x2)" (1) 


This foi mula is repeated on Fig. 2. 


Figure 3 shows three sample keel sets that were used to estimate scattering. They will 
be identified in the -- mainder of this report as sets A, B, and C. Sigma indicated or. each one is 
the standard deviation of ice depth in meters. As can be seen, there is considerable variation in 
this parameter, although the three sets were generated by selecting random values from the 
same distrioutions. If the three sets are combined into one long set, it has a sigma of 2.1 m, or 
nearly average tor observed Arctic ice roughness. 
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Figure 3. Three keel sets generated by random numbers. The 
standard deviation of ice depth is given on each set. 
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Ray scattering functions were derived as follows. A number of rays were launched at 
tne same angle, but at successive ranges to simulate a plane wave directed at the keel set ata 
given angle of incidence. Usually 800 rays were used. Rays leaving the keel set and crossing a 
given receiver depth are sorted into angular bias. Any reflection of transmission losses at 
boundanes aiters the intensities of the accumulated rays. Spreading loss is accounted for by the 
collection into angular bins. Rays are equally spaced along the plane wave front and therefore 
represent equal energy. The sum cf cay contributions in any bin is a measure of relative 
intensity. 


Two scattering processes are allowed in the computations: reflection at the ice-water 
interface anc refraction through the keel. Refraction is permitted at both the shear and 
compressional speed 17 the ice. However, few rays strike the keels at sufficiently large grazing 
angles for transmission at the compressional speed, so shear wave paths through the keels 
predominate. Edge-diffracted rays were included in Ref. |. But since their contnbution to the 
field is small, they are not included here. 


The computer program used here requires a non-zero sound speed gradient, so the rays 
will have curvature. Since our purpose here is to determine the scattering properties of the set 
of keels, we prefer straight line rays. The curvature is m:nimized by using a smal! gradient and 
by placing source and receiver depths just below the keels. These depths were usually !0 or 
20 m. depending on the deepest keels in the set. A sound-speed gradient ot 0.0002 sec! was 
used in all computations here. 


The retlection and transmission properties of ice can be computed following Ewing, 
Jardetsky, and Press (Ref. 4). The acoustic parameters required are the sound speed in the 
water, the shear and compressional sound speeds in the ice, and the density of water and ice. 
Table | shows the values of these parameters used here. These are typical observed values for 
Arctic sea ice. The shear speed varies over a large range of values, depending on the condition 
of the ice. Observed values easily vary from 1400 to 2200 m/s. 


Table |. Typical acoustic parameters in Arctic sea ice and water. 


Parameter Value 
Sond speed in water 1440 m/s 
Compressicnal sound speed in ice 3600 m/s 
Shear cund speed in ice 1600 m/s 
Water density 1026 kg m3 
Ice density 926 kg) m3 


SCATTERING RESULTS 


Figure 4 shows the angular distribution of scattered rays. The sum of rays !n each 
2-degree angular bin '!s expressed as a decibel by taking -10 log of ‘he ratio of the ray sum ina 
bin to the total! rays launched. Keel distmbutior A of Fig. 3 was used. The svrprising result here 
is that almost all rays are scattered into the forvars direction. Very small backscattering lobes 
can be seen for the steeper rays. Note also that a significant specular component does not show 
up for angles of incidence below 6 degrees. 
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Figure 4. Energy scattered into 2-degree bins as a fraction of 
total energy in a piane wave striking keel set A at the designated 
angle (grazing ang'e for a ray). 
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T12 smal. backscatter is an urrea’istic result because of the limited slope angles of the 
parabolic keels “his wil! be discussed .11 later sec‘1ns, which give detailed analyses of the 
¢ Mmehanic: crea inz these curves. 11 generia:, the ft:dominance of forward scatter over 
dackscatter re‘.its frorm the wafer-ice r:flection : oefficients. One set of plane wave reflection 
and trans niss C4 coet,icien's is sFown in Fig. 5. There are similar sets of reflection coefficients 
for the shear and compressional waves impinging cn the ice-water interface. The sea ice shear 
sourd spezd of 16C0 in/s from Tatle | gives a critival grazing, angie of 26 degrees. Below this 
angle, pervect reflection ccew’s for iossless ined.a. Above this angle, the reflection loss increases 
rapid’y. Simettare.: islv the * unsm ss on Icss into sear waves decreases rapidly. The ice 
becorces very triaas*arent tc ihe sounc. Ths t-anspz rency at high grazing angles is a major 
cause of the prefe-e‘itial forward scatter of the rays. 
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Figure 5. Piane wave reflection and transmission coefficients at 
@ water-ice interface for a compressiona! wave in the water. 


Figure 6 shows angular scattering diagrams for three angles for the three keel sets of 
Fig. 3. Set A has a considerably higher roughness factor, or sigma, than sets B or C. However, 
no distinctly different characteristics are apparent in set A. 


* 
al 


3 Figure 7 gives an idea of what art of the scattered energy is contributed by reflection 
oi and what par by refraction through the keels. Figure 7 is equivalent to panel A in Fig. 6, but 
Rs only the reflected rays are included. As can be seen, the refracted energy accounts for a 
j oa) predominant part of the scattered field. This refracted energy tends to remain closer to the 
i se specular direction, with the reflected energy accounting for much of the higher angle scatter. 
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The specular reflection coefficient is the product of this investigation that is most 
immediately usable in propagation loss programs. It essentially measures the energy that 
remains coherent upon reflection from the rough surface. Figure 8 shows this reflection 
coefficient as a function of grazing angle for the three keel sets. Below a grazing cngle of 6 
degrees, the larger roughness of set A seems to produce a larger scattering loss. At larger 


x 


ns re 5 re ree ee Se ere er ee 


angles, the roughness seems immaterial. The rays that contnbute to this specular reflection are 
almost all those that reflect from the smooth surface without touching a keel. The function 
shown here is thus a measure of the unshadowed surface as seen from both the source and 
receiver. 


Theory tor constant-depth bosses holds that the reflection loss, which becomes large at 
_ low grazing angles, should recover at very low angles and become zero at zero-degree grazing 
angles. However, computations with this program down to 0.75 degree have shown no 
tendency towards recovery of the reflection loss. ; 
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Figure 6. Angular scattering functions by 2-degree bins for keel 
sets A, B, and C. 
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Figure 7. Angular scattering function 
for keel set A, where only reflected rays 
are used. Transmitted rays are omitted. 
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Figure 8. Specular scattering coefficients as a function of 
grazing angle for the three keel sets These are equivalent to the 
pvaticcing function values in the 2-degree bin containing the 
specular angie. 
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GEOMETRIC IMPLICATIONS OF SCATTERING ANGLES 


In this section the general shape of the scattering functions will be discussed. In partic- 
ular, the backscattering function will be shown to be underestimated because of the limited 
slopes of the ice model. This topic is discussed first because it is important that the backscatter 
functions of this report not be used for estimation or modeling purposes. A geometric discus- 
sion of the specular reflection function and the scattering curves will follow to explain some of 
the features of these curves. 


BACKSCATTER 


Two notable features of such scattering curves as those of Fig. 4 and 6 are the null in 
the downward direction and the small backscattering lobes. Both these features result from an 
interplay between the parabolic shape and the reflection coefficient. Figure 5 shows reflection 
and transmission losses for a water-:ce interface. These are computed by the method of Ewing, 
Jardetsky, and Press (Ref. 4). Table | gives the acoustic parameters used. Computations of the 
reflection loss for a wide range of these parameters are given by Mayer, Behravesh, and Plona 
(Ref. 5). Although the losses in Fig. 5 are plotted as functions of the phase velocity, an 
equivalent grazing angle scale is given at the top of the figure. 


As a result of the acoustic parameters chosen here, acoustic energy begins to penetrate 
the ice-water interface at a grazing angle of about 26 degrees. The next increase in grazing 
angle of a few degrees finds almost ail energy transmitted into the ice via the shear wave and 
almost none reflecting. At an angle of 45 degrees, less than 1% of the energy is reflected. A ray 
reflecting from a sloping surface is diverted by twice the slope of the surface. (Exact expres- 
sions will be given later.) The rays scattered downward near 90 degrees have reflected from a 
slope near 45 degrees and hence suffered a large reflection loss. The nulls in the downward 
direction on the scattering curves are images of the null in the reflection coefficient of Fig. 5. 


Near a grazing angle of 66 degrees the compressional wave in the ice is excited. From 
this angle on to 90 degrees 10% or more of the energy is reflected. Figure 5 shows good 
reflection coefficients starting at about 60 degrees. Such slopes will reflect at 120 to 180 degrees 
or backscatter the energy. However, the parabolic ice keels used in this study have maximum 
slopes of 51.3 degrees. Therefore only a few of the rays of highest incident angles can reach 
these good reflectors for backscatter. 


To test the assumption that higher slope angles would produce more backscatter, a run 
was made with keels of half width. This was a relatively simple test to make. The keels of set A 
were reduced to half their width, with all other parameters remaining the same. This increases 
the maximum slope angle from 51.3 to 68.2 degrees. This should permit many more rays to fall 
into the preferred grazing-angle range. The additional 17 degrees should increase the backscat- 
tered angular bins to angles about 34 degrees higher. Figure 9 shows the results, leaving little 
doubt that the higher slope angles do control the backscatter and that if all slopes up to 90 de- 
grees occurred in the model, the backscattered field would be continued completely around to 
180 degrees. The accuracy of such a field would depend on the overall accuracy of the ice keel 
model. 


Figure 9 shows that slopes above 51.3 degrees increase the backscattered field. How- 
ever, the forward-scattered fields should not be accepted in the place of those of Fig. 6 because 
the modeled keels do not have dimensions equal to the best estimates of real ice keels. 
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Figure 9. Angular scattering functions 
for a set similar to keel set A, except 
that the keel width is reduced by half 


10 deg 


SURFACE SHADOW!NG 


The curves of Fig. 8 appear to be a measure of the unshadowed surface. In the 
following Ciscussion, a theoretical derivation of the shadowing under simplifying assumptions 
will be derived and compared with the model results. Let us remove the randomness and 
assume that the keels are evenly distributed, 9.5/km, and have a depth of 4 m. These are the 
means of the uniform and Rayleigh distributions from which the random sets were made. We 
need to compute for plane wave angle 6 the fraction of surface not shadowed by the keels, f, 
and convert to decibels by taking -10 log f& A general expression will be given first. Figure 10 
shows the wave front, s, for a ray angle of 30 degrees for a specific keel. Rays that reflect in 
front of and then strike the keel are represented by use of a reflection of the keel on the surface. 


Let d be the depth. b the half width of the keel, and s; the length along the approaching 
piane wave front that will intersect the keel. Then 


s;=[ (b2 tan? @ + 46d tan 6 + 4d2) cos @)/4d (2) 


The width of wave front that intersects the keel after reflection in front of the keel is given by 
changing the sign of 6. 


Thus the total length of wave front intercepted, s,, is 
5, =( (b2 tan2 6 + 4a2)cos 6] / 2d (3) 


The surface length shadowed by this length of wave front is s,/ sin 6. These equations 
apply up to the maximum slope angle of the keel, max. OF 


max = arctan(2d/ b) (4) 


At greater values of 6, the shadow is just the width of the k:el. 
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Figure 10. Portion of a plane wave front intersecting a keel 
and its reflection on the surface. A coordinate axis labe!ed s 
is centered on this wave front. 


The half width 6 was 1.6 d in the original distribution. This gives a maximum slope 
angle of 51.34 degrees. Since the keels were randomly wiuened by Eq. ! to simulate random 
directional orientation of the keels, a larger value of b is appropriate. The mean widening 
factor can be determined by integrating Eq. | from 0 to |. A value of 1/2 is obtained. This 
gives an average half width of 10.05 m, compared ith 6.40 m without widening. The percent 
of shadow for b = 10.05 m and converted to decibels as described above is plotted in Fig. 11. 
The ray theory values from Fig. 8 are included for comparison. The shadow covers the 
complete surface at a grazing angle of 4.39 degrees, and the loss is infinite. At alrnost all angles 
above 8 degrees this surface shadow curve lies within the ray theory data points. The smaller 
loss of the ray theory over the shadowing curve for angies below 7 degrees may result partly 
from scattered rays falling into the 2-degree bins in which rays are collected. However, the 
largest difference arises from the randomness of che keels, which do not shadow the surface as 
efficiently as do keels of uniform size and spacing. 


A quick correction for this =ffect is to assurre that the shadows are randomly spaced 
and subtract the fraction of double-shadowing from the fraction of shadowed surface, giving 
x2/2, where x is the probability of a point on the surface being shadowed. This shadowed 
fraction is in turn corrected for triple shadowing, <.d so on. This series is a representation for 
the exponential function, a familiar result from prod bility theory. Thus the unshadowed 
surface |-x is represented by exp (-x). The lower line in Fig. 11 shows this loss for a randomly 
shadowed surface. The resulting curve seems to serve as a lower limit to the reflection loss. 


Had we used the narrower keels having a half width of 6.4 m. Se upper line would 
have been very similay. It is about 0.22 dB less at 5 degrees, 0.10 dB less at '? degrees, and 
back up to 0.15 dB less at 20 degrees. These narrower keels fully shadow the su::2ce at a 
4.37-degree grazing angle. The shadowing strength of the keels depends very little on their 
width, but rather on depth and spacing. 
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Figure 11. Estimated specular reflection coefficients are 

compared with those determined by accumulating rays. The 
upper estimate is based on evenly spaced, similar keels. The 
bottom estimate is based on randomly spaced, similar keels. 


SCATTERING FUNCTIONS 


The scattering functions are a sum of the reflected and refracted fields. A full analytical 
description of neither of these fields is attempted in the following discussion. However, some 
equations and plots are given to clarify some aspect of the fields. 


Returning to Fig. 10, the slope of the parabola at the point a ray strikes can be 
determined. The ray is defined by the value of s or distance from the center of the wavefront 
that intersects the keel. The angle of the wave front, 6, is defined as the angle the rays make to 
the horizontal, where the angles of rays arriving from below the horizontal are positive. The 
slope of the parabola ¢ is given by 


$ = tan-![ -2ar/(b2)] (5) 
where the intercept range, r, is given by 

r= { (tan @)/(d/b=) -[ (tan2 6)/(d2/b4) + 4d - s/cos @)/(d/b)'2}/ 2 

Figure 12 summarizes the use of these last two equations for the keel of Fig. 10. It isa 
plot of certain slope values at points of intersection between a ray and the keel. The rays are 
determined by ray angle 6 and s, the position on the wave front as defined in Fig. 10. The two 


lines labeled “limiting ray” are the tangent rays, and they intersect the parabola where its slope 
is -8. The angle between ray and parabola is 6 + @ or zero at the point of tangency. The width 
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Figure 12. Locus of points along the wave front from which rays 
of a given grazing angle strike the keel at designated angles. 


between these limiting rays increases by only 0.4 m from a ray angie of 0 to 30 degrees. The 
maximum slope of 51.3 degrees occurs where the keel meets the surface, so the ray reaching 
this point is the boundary between the direct rays and the surface reflected rays. As can be seen 
from the 51.3-degree line in Fig. 12, as the grazing angle increases, a growing proportion of the 
wave front intersects the keel directly and a decreasing portion intersects the keel after 
reflection from the surface. The proportion would reach zero if the figure were extended to a 
ray angle of 51.3 degrees. 


The slope of 26 degrees is of interest because, as was discussed eaclier, at grazing angles 
of 26 to 60 degrees the rays are largely transmitted into the ice. The locations of rays that 
intersect the keel where its slope, or @, is 26 degrees are labeled 26-degree slopes on the figure. 
However, of more interest are the rays for which 6 + @ is 26 degrees. This is the true dividing 
line between reflected and transmitted rays. Rays with this property are shown by dashed lines 
in the figure. Note that only about 1$% of the rays are perfectly reflected. Most of the rest meet 
the keel at angles that favor transmission through the keel a. a shear wave. However, at an 
incident angle above 60 degrees there is a more even split between reflected and transmitted 
rays. The wedge between the lines labeied “$1.3-deg slope” and “60-deg incident angle” repre- 
sents these rays, which comprise more than half the wave front at ray angles above 25 degrees. 
The reflected energy from these rays will be directed backward to the backscattered region. 

The direction that a reflected ray leaves the keel, 9,, is given by @ + 2d. Figure 13 shows 
these reflected angles for an incident wave front of 20 degrees. The line labeled “slope angle” is 
the slope on the keel and is equivalent to the angles along the 20-degree abscissa in Fig. 12. The 
reflected angle is then twice this slope plus 20 degrees for the direct rays on the nght side of the 
figure. At a wave front location of ~2.2 m, the rays become surface-reflected, and 6 is 
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Figure 13. Angles at which rays emerge after striking or passing 
through the keel of Fig. 10, as a function of position on the wave 
front. The slope of the keel at the point of contact is labeled 
“slope angle.” 


-20 degrees. This causes the vertical discontinuity of 40 degrees in the reflected angle. Reflected 
angles less than 0 degrees are headed upward after reflection and reflect from the surface 
behind the keel, becoming downgoing or positive angles of the same size. 


The final angle at which refra_ted rays emerge from ti ¢ keel cannot be easily 
characterized. Figure 14 shows the change in angle at the intertace for the water wave to ice 
sear wave interaction for the sound speeds of Table |. A 45-degree line is drawn in to indicate 
the size of the angle change. At high grazing angles or near-normal incidence, very little change 
occurs. These angular changes occur at transmission both into the keel and from it. The 
additiona! change in ray angle is a result of change in slope of the keel at the entering and 
exiting points. The ray may also reflect from the flat upper surface of the keel model. 


A compressional wave begins when the grazing angle equals 66 degrees and ends when 
the rays reach the upper corner of the keel. The upper end of this line has a small transmission 
coefficient, hut towards the lower end provides occasional downward-directed rays, as can be 
seen in Fig. 4 and other such scatter diagrams. 


The accumulation of the above effects produces exiting angles for a 20-degree ray sei, 
shown by the dashed lines in Fig. 13. At the critical grazing angle near 26 degrees, the reflected 
and refracted rays correspond in angle. The dashed line therefore meets the reflected angle 
curve at the points where the slope angle is 26 + 20 degrees. The nearly flat section through the 
middle of the refracted angle line is from rays that reflect from the upper surface of the keel. 
The discontinuity at either end cf this section is the ray that intersects either the leading or 
trailing corner of the keel and surface. 
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Figure 14. Refraction angle at water-ice interface for the water 
wave and ice shear wave. 


i 
a 
J The final scattering function for a single keel is the convolution of the angles of Fig. 13 
with the reflection losses of Fig. 5. However, the transmission losses of Fig. 5 contain ray- 

fe} spreading terms that must be removed, since equal intervals along the s scale of Fig. 12 
represent equal energy, and spreading effects are built in. 


‘ay 

¥ MODEL SHORTCOMINGS 

Me . The models used here fail to represent reality in many ways. Some of their 

ity shortcomings are listed here. The results of limitations on slope angles were discussed in the 
sty : section on backscatter. Perhaps the next most noticeable shortcoming is the absence of 

Bs reflection from the far sides of the keels. These reflections might add significantly to the 

a backscattering strengths of the keels. In a pure ray construction treatment, as used here, one 
ef must terminate internal reflections at some number or level. Here, they were terminated at 
te zero. 

aft A third important shortcoming is the absence of an ice plate which could alter ray 


paths. This shortcoming does not cause serious angular dislocations though, because the 
transmission coefficient is such that rays directed into the ice plate will readily pass back into 
oe the water. A nonflat upper ice surface, however, will disturb ray paths. 
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Finally, attenuation of shear waves in the ice is omitted. This may very well add a 
decibel or more loss to paths through the ice, changing the angular scattering diagram. This 
effect and others upon the reflection coefficient at the ice plate might reduce the specular 
reflection coefficient, which depends upon reflection from the plate. 


CONCLUSIONS 


A collection of ice keels has been modeled by parabolas that have distributions of size 
and spacing that match estimates of real keels. Rays representing a plane wave have been scat- 
tered from this set of keels, with resultant scattering angles and specular reflection coefficients. 


The reflection-transmission charactenistics of the water-ice interface are shown to be 
important in determining the angles into which energy is scattered. The particular geometric 
properties of the parabolas also have a strong effect on the estimated scattering angles. This is 
particularly true of backscatter, which is underestimated by the parabolic shapes. 


Specular reflection coefficients from the keels are estimated by computing the ratio of 
shadowed surface io total surface. 


For grazing angles up to 20 degrees. transmission through the keels accounts for more 
scattered energy than reflection trom the keel. Below 20 degrees, this transmission is 
predominantly by shear waves rather than compressional waves. 


RECOMMENDATIONS 


The ray theory program utilized here permits keel shapes to be approximated by 
several parabolic or straight-line segments, although only a single parabolic segment per keel 
was used here. Such compound shapes could more nearly model real keels. Pecause the shapes 
used here lead to defective estimates of backscatter, the use of compound shapes is particularly 
indicated for keel backscatter estimations. 


Before more appropmate keel shapes are generated, much information is needed about 
detailed keel shape and structure. If keels tend to remain collections of independent ice blocks 
or slabs throughout their life, then the current analysis, which is based on solid ice shapes, has 
very limited validity. If keels do solidify and if their shape is modified by local melting and 
freezing, then a solid ice model can be appropriate. Clearly, a great deal of observed 
information on shapes and slopes is needed to produce random model shapes that match 
readiy. Upsounder echo traces might serve as sources of this needed information. 


If acoustic paths through !ce keels are important, as the current analysis indicates, then 
the internal elastic-acoustic properties and degree of homogeneity of keels is irmportant. Of 
particular importance are the sound speed and the attenuation of the two wave types, shear 
and compressional. A systematic study of the ice in keels is needed. 


Acoustic experiments performed in the vicinity of a well-observed keel could reveal 
much about the acoustic properties of keels. The basic requirements are receivers at several 
depths, both in front of and behind the keel, and a vanable source angle. The ability to 
separate multipath arnvals would also be useful. 
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